Generalized van der Waals theory of liquid-liquid phase transitions 
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In the framework of the thermodynamic perturbation theory for fluids we study how the phase 
diagram of an isotropic repulsive soft-core attractive potential, where a liquid-liquid phase transition 
exists in addition to the standard gas-liquid phase transition, changes by varying the parameters of 
the potential. We show that existence of the liquid-liquid transition is determined by the interplay 
of the parameters of the potential and the structure of a reference liquid. 
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Despite the growing interest to the possible polymor- 
phic phase transitions in liquids and glasses pi Q the 
nature of different phases which can be found in dense 
(and possibly metastable) liquids is still puzzling. The 
coexistence of different phases, while common for mix- 
tures, is unexpected for a simple fluid. In principle, the 
rules of thermodynamics do not forbid the existence of 
more than two different fluid phases in a simple fluid, 
however, from a common experience only two of them 
are well known: a low density fluid (vapor) and a high 
density fluid (liquid) phases. At the same time in recent 
years experimental evidences of such features of phase 
diagram as liqiud-liquid transitions, polyamorphism, etc 
appeared for a wide range of systems including water, 
Si,I,Se,S,P, etc 0, S Ull SBS S H The 
complexity of the phase diagrams in these substances 
may be a result of complex interactions depending on 
the intermolecular orientations. 

At the same time exploring the possibility that simple 
fluids interacting through isotropic potentials may ex- 
hibit the similar behavior represents a serious challenge 
for theorists. 

It was shown recently through molecular dynamics 
simulations that a system of particles with the isotropic 
repulsive soft-core attractive potential may have high- 
density and low-density liquid phases |2(| ■ This po- 
tential may be considered as an effective potential result- 
ing from an average over the angular degrees of freedom 
for systems where the position of the minimum approach 
between particles depends on their relative orientations 
like in the case of the hydrogen bonding between water 
molecules Q, H, 0, |23| . This potential may be also used 
to model interactions in a variety of systems including 
liquid metals, colloids, silica [lll23|. 

After the pioneering work by Hemmer and Stell |16| . 
where the soft core potential with an attractive interac- 
tion at large distances was first proposed for the quali- 
tative explanation of the solid-solid critical point in ma- 
terials such as Ce or Cs, a lot of attention was paid to 
the investigation of the properties of the systems with 
the potentials that have a region of negative curvature 
in their repulsive core. In spite of the simplicity of the 
model, the physical mechanism that causes the liquid- 
liquid phase transition in such systems is not completely 
understood. As was emphasized in [24l |25| it arises from 
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FIG. 1: An isotropic potential with an attractive well and 
a repulsive component consisting of a hard core plus a finite 
shoulder. 



an interplay of the different components of the pair in- 
teraction. In Refs. pil l2o| authors tried to disentangle 
the role of each component to investigate the dependence 
of the phase diagram on the potential parameters. In 
Ref. [24J the results of molecular dynamics calculations 
performed for several sets of parameters were presented. 
The resulting behavior of the critical points was inter- 
preted through a modified van der Waals equation where 
the effect of the repulsive shoulder at different densities 
p and temperatures T can be taken into account by an 
effective excluded volume depending on both p and T. 
In Ref. [25| the same analysis was undertaken by using 
an integral equation approach in the hypernetted-chain 
approximation. In Refs. 0, HH ^ was shown that the 
high-density liquid-liquid critical point can be found only 
when there is some kind of balance between the attractive 
and repulsive parts of the potential. 

It should be noted that it is widely believed (see, for 
example, [2H W\ ) that a fluid-fluid transition should be 
related to the attractive part of the potential. However, 
in Ref. psj in the framework of the generalized van der 
Waals theory it was shown that the purely repulsive step 
potential is sufficient to explain the high density liquid- 
liquid phase transition. In the present paper we apply 
this theory to the study of the phase diagram of the sys- 



2 



T=5.0 

solid o=1.2 
dash a=1.4 
dot o=1.6 
dash-dot o=1.8 



0,3 0,4 o,s 0,6 

n 



FIG. 2: The pressure isotherms for the system with e/ei = 5 
for different values of a/d for T = 5 (ai = 2). 



tem of particles interacting through an isotropic potential 
with an attractive well and a repulsive component con- 
sisting of a hard core plus a finite shoulder. This potential 
can be written in the form: 



$(r) 



00, r < d 

Ei, d < r < a 

£2, a < r < o\ 

0, r > ai 



(1) 



The potential is shown in Fig. QJ. We apply to the 
problem the first order thermodynamic perturbation the- 
ory for fluids. The soft core of the potential is treated 
as perturbation with respect to the hard sphere poten- 
tial. In this case the free energy of the system may be 
written in the form [2^, |2jJ : 



F 



f hs _ 1 a 
Nk B T 2 Pf 



ui(r)g H s(r)dr, 



(2) 



where p = V/N is the mean number density, = 1/kgT, 
u\{r) is the perturbation part of the potential u\(r) = 
$(r) — $Hs(r), <&Hs{ r ) is the hard sphere singular po- 
tential, gHs( r ) is the hard sphere radial distribution func- 
tion, which is taken in the Percus-Yevick approximation 

m 

To calculate Fhs, one can use, for example, the ap- 
proximate equation j3fj: 



Fhs „, . . , , , 477 - 3T7 2 
^__ = 31nA-l + lnp +Tr — _. 



(3) 



Here A = h/{2-Kmk B T) 1 / 2 and 77 = 7rpa 3 /6. 

Further in this paper we use the dimensionless quan- 
tities: f = r/d, P = Pa 3 /e, V = V/Na 3 = 1/p, 
T = ksT/e, omitting the tilde marks. 
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FIG. 3: Phase diagram of the system of particles interacting 
through the potential (1) for different values of a jd and fixed 
value of ai/d = 3. Solid lines correspond to binodals and the 
dashed lines - to spinodals. Inserts show the characteristic 
isotherms. 



Results of calculations are demonstrated in Figs. 121 II 
In Fig. [2] a family of pressure isotherms is shown for the 
system with e/e\ = 5 for different values of a/d for T = 5. 
The van der Waals loops in the isotherms at low temper- 
atures are clearly seen, this indicates the existence of the 
first order liquid-gas and liquid-liquid phase transitions. 

Using the Maxwell construction we are able to calcu- 
late the equilibrium lines of the liquid-liquid phase tran- 
sitions at different values of a/d (Figs, l-il 11 (binodals) 
which are shown as solid lines in the figures. The dashed 
lines correspond to spinodals which as usual are calcu- 
lated from the condition dP/dp = 0. In Fig. [3] the 
evolution of the phase diagram is shown as a function of 
a/d for the fixed value of a\/d = 3. We can see that the 
phase diagram starts from one gas-liquid transition for 
a/d = 1.3. If value of a/d increases the second (liquid- 
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FIG. 4: Phase diagram of the system of particles interacting 
through the potential (1) for different values of a/d and fixed 
value of Aa = a±/d — a/d — 0.5. Solid lines correspond to 
binodals and the dashed lines - to spinodals. Inserts show the 
characteristic isotherms. 



liquid) transition develops. When difference ai/d — a/d 
is small enough (« 0.2) the gas-liquid transition disap- 
pears, and one has the only liquid-liquid transition, as 
was discussed in Ref. j2^|. We cannot extend the tran- 
sition lines down to zero temperature and make calcula- 
tions for u\/d < 1.3 because of limitation of the pertur- 
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FIG. 5: The potential @J for r = 1.1 (line 1) and r = 1.16 
(line 2). The line 3 is obtained from the line 2 by the shift 
downward on the height of the first maximum. 



bation approach. 

Fig. 01 shows the evolution of the phase diagram of the 
system of particles interacting through the potential (1) 
for different values of a/d and fixed value of the differ- 
ence Act = a\/ d — a / d = 0.5. It should be noted that 
in contrast to Fig. |31 the liquid-liquid phase transition 
disappear for values of a in the vicinity of a « 2. This 
seems to contradict the assertion in Refs. p3. 125} that 
the existence of a liquid-liquid phase transition is simply 
determined by the some kind of a balance between the 
repulsive and attractive parts of the potential, but de- 
pends also on the positions of the maxima of a reference 
radial distribution function (or the structure of the refer- 
ence liquid). This is seen from the right hand side term 
ofEq. 0. 

To illustrate the application of the generalized van der 
Waals theory to other systems let us consider the poten- 
tial of the form: 



aexp(— ar) cos(2fc/(r — Tq)) 



+ b 



18 



(4) 



where a = 0.1, a = 0.5, kf — 4.1, a = 0.331, b = 
0.42 x 10 s . This potential may be used, for example, 
for a qualitative modelling of effective potentials of some 
metals j33|. In Fig. [S] this potential is shown for ro = 1.1 
(line 1) and r = 1.16 (line 2). The line 3 is obtained 
from the line 2 by the shift downward on the height of 
the first maximum. 

In Fig. the corresponding families of isotherms are 
shown for different temperatures. To calculate the equa- 
tion of state corresponding to the potential (QJ we ap- 
ply the thermodynamic perturbation theory. The second 
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part of the potential was considered as a reference 
system. 

*°M = ^) 18 - (5) 

Potential §o(r) is approximated by the hard sphere po- 
tential with an effective diameter which depends on the 
density and the temperature [33|. 

From Fig. [5] one can see that there are only high- 
density liquid-liquid transitions for the potentials corre- 
sponding to the lines (l)and (2) in Fig. (upper fig- 
ure). These transitions exist due to negative curvature 
of the potential Q. There are no gas-liquid transitions 
in these cases because the strength of an attraction is 
insufficient. For the potential corresponding to the line 
(3) in Fig. [5] there is also a low-temperature gas-liquid 
transition shown in the insert in the lower figure. 

We would like to emphasize that we do not claim that 
the first order perturbation scheme, which was used in 
the present article, gives the high precision quantitative 
results, however, it seems reliable enough to give correct 
qualitative description of the liquid-liquid transition in 
the system with the potentials Q and (0J. It should be 
noted that the second-order perturbation theory gives 
qualitatively the same results. 

Finally, let us make some remarks on another possible 
mechanism of a liquid-liquid phase transition. The crys- 
talline solid state is characterized by a long-range posi- 
tional order of the atomic density and associated bond 
orientational order belonging to one of the well-known 
allowed crystalline symmetries. In a diffraction exper- 
iment, S- function Bragg reflections are observed in the 
structure factor 5(q) that shows the associated point 
symmetry of the lattice. In the case of the long range 
bond-orientational order, one can observe the modulation 
of the structure factor which corresponds to the symme- 
try of the nearest neighbor environment of a particle (see, 
for example, [34j]). 

However, one can imagine the situation when the 
isotropy of three or four-particle correlation function is 
broken, but the symmetry of one and two-particle corre- 
lation functions is unchanged. In this case the structure 
factor of the system is the same as for the isotropic liquid 
because it depends only on the two-particle correlation 
function. This is why we call this type of the symmetry a 
" hidden symmetry" - it can not be detected in diffraction 
experiments. 

So, using the diffraction experiments it is impossible 
to distinguish a liquid-liquid transition without changing 
the symmetry of the correlation functions and the transi- 
tion accompanied by the breaking of a hidden symmetry. 
It seems that the only way to answer the question what 
kind of transition one observes in an experiment is to 
analyze the whole phase diagram. In principle the hid- 
den symmetry transition may be of the first order, but in 
contrast to the liquid-liquid transition discussed in this 
article, it separates the phases with different symmetries 
and can not end at the critical point. There are two pos- 
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FIG. 6: The families of isotherms corresponding to the po- 
tentials in Fig. There are only liquid-liquid transitions for 
the potentials corresponding to the lines (l)and (2) in Fig. 
(upper figure). For the potential corresponding to the line (3) 
in Fig.^Jthere is also a low-temperature gas-liquid transition 
shown in the insert in the lower figure. 

sibilities for the hidden symmetry phase transitional) 
it may end at the tricritical point, and measurements 
above the tricritical point should reveal the thermody- 
namic anomalies but no density change; (2) it may inter- 
sect the gas-liquid transition line. In this case one can 
expect a sharp bend on the gas-liquid transition line at 
the point of intersection. 

To our knowledge there are no calculations exploring 
the idea of breaking of the symmetry of higher order 
distribution functions in application to the problem of 
liquid-liquid transitions. 
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